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Abstract 

We present if -band polarimetric imagery of UX Tau A taken with HiCIAO / AO 188 
on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk 
object, with a gap structure separating its inner and outer disks. Our imagery taken 
with the 0.15" (21 AU) radius coronagraphic mask has revealed a strongly polarized 
circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial 
resolution of 0.1" (14 AU). It is inclined by 46° ± 2° as the west side is nearest. 



Although SED modeling and sub-millimeter imagery suggested the presence of a gap 
in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, 
we detect no evidence of a gap at the limit of our inner working angle (23AU) at the 
near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) 
to light scattering by dust grains in the disk. However, neither polarization models of 
the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations 
were consistent with the observed azimuthal profile of the polarization degrees of the 
disk. Instead, a geometric optics model of the disk with nonspherical grains with the 
radii of 30 /im is consistent with the observed profile. We suggest that the dust grains 
have experienced frequent collisional coagulations and have grown in the circumstellar 
disk of UX Tau A. 

Key words: (stars:) planetary systems: protoplanetary disk - techniques: high 
angular resolution - techniques: polarimetric 

1. Introduction 

A protoplanetary disk around a young stellar object is the site of planetary formation. 
The core accretion model (e.g., Nakagawa et al. 1983) and the gravitational instability model 
(e.g., Boss 1998) are two possible methods by which the process of planet formation might 
occur. In the core accretion model, a key process is dust grain growth by collisional coagulations. 
Several fundamental processes of the dust growth have been proposed under simple assumptions. 
Dust grains with sizes less than a few tens of micrometers orbit a central star with the same 
velocity as gas (Adachi et al. 1976). They frequently collide and coagulate with each other. 
Simultaneously they begin to settle into an equatorial plane of the disk (Nakagawa et al. 1981). 
Dust grains with different settling velocities also experience collisional coagulations. As a result, 
planetesimals (radii ~ 10 - 100 km; Kokubo & Ida 1998) are formed in the mid-plane of the 
disk. Bodies between a few tens of micrometers and a few kilometers are dominated by gas 
drag. Such bodies are thought to migrate towards the central star in a short timescale (e.g., 
Nakagawa et al. 1986). This rapid inward migration poses a challenge to planet formation 
theory. By contrast, for planetesimals larger than a few kilometers, gravitational interactions 
are significant and relative velocities between the bodies increase due to mutual perturbations. 
This situation allows faster growth of the larger bodies and leads to formation of planetary 
embryos (radii ~ 1000 km; Kokubo & Ida 2000). The core accretion model is one of the 
plausible model of planetary formation. However it is still debated how sub-micron sized grains 
are transformed into kilometer sized bodies. 

An infrared excess in the spectral energy distribution of a young stellar object provides 
indirect evidence for the presence of a circumstellar disk. The spectral energy distributions 
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of T Tauri stars with continuous, optically thick disks have strong infrared excesses in the 
near-infrared to far-infrared wavelengths (Williams & Cieza 2011). Recently, objects with large 
continuum excesses in the mid- to far-infrared wavelengths but no excess in the near-infrared 
wavelengths have been discovered, called 'transitional disk objects' (TDOs; Calvet et al. 2002, 
Calvet et al. 2005). Such an object is considered to lie in the transition state between classical 
T Tauri stars (CTTSs) and weak-line T Tauri stars (WTTSs). It is expected that the inside of 
the disk has been cleared out by dust accumulations and/or formations of protoplanets. Photo- 
evaporation process is also proposed as a mechanism to create a transitional disk. On the other 
hand, 'pre-transitional disk objects' (PTDOs) are considered to be in the evolutionary phase 
before reaching the transitional disk phase. They have mid- to far-infrared excesses similar 
to TDOs, however, they also show small excess in the near-infrared wavelengths (Espaillat 
et al. 2007). The slight near-infrared excess implies that optically thick material remains in 
the innermost part of the disk. Espaillat et al. (2010) proposed that the planetary formation 
process is expected to progress in the gap structures between the inner and outer disks. 

High spatial resolution coronagraphic imaging polarimetry is one way to directly diagnose 
young circumstellar disks. The degree of polarization depends on scattering angle, grain size, 
and composition. By constructing the polarization profile of the disk, we are able to investigate 
size and composition of dust grains. Silber et al. (2000) conducted infrared polarization imaging 
observations of a circumbinary disk around GG Tau with NICMOS mounted on the Hubble 
Space Telescope. The circumbinary disk shows strong polarization degrees of ~ 50% at 1 /zm 
wavelength. The polarization azimuthal profile indicated Rayleigh-like scattering from dust 
grains with sub- micron size. Hashimoto et al. (2011) conducted if-band polarization imaging 
observations of AB Aur with HiCIAO/A0188 on the Subaru Telescope and revealed spiral 
structure in the outer part and the double ring structure at the inner part of the circumstellar 
disk. There was a number of studies in the last years, using the same observational technique 
for similar sources, e.g., Apai et al. (2004; TW Hya), Oppenheimer et al. (2008; AB Aur), 
Perrin et al. (2009; AB Aur), Quanz et al. (2011; HD100546), Quanz et al. (2012; HD97048). 

UX Tau is a T Tauri multiple system (Jones & Herbig 1979) in the Taurus molecular 
cloud (distance ~ 140 pc; Elias 1978). It consists of a primary star (UX Tau A), with UX 
Tau B separated by 5.86" and UX Tau C separated by 2.63" from the primary. UX Tau B 
is itself a binary system with the separation of 0.14". The spectral type of UX Tau A is K2 
(Kraus & Hillenbrand 2009). Its spectral energy distribution shows a slight excess in the near- 
infrared wavelengths and significant excesses in the mid- and far-infrared wavelengths. These 
characteristics indicate that UX Tau A has an optically thick inner disk separated from an 
optically thick outer disk by a gap, i.e. it is a pre-transitional disk object (Espaillat et al. 
2010). Model fits of the SED suggested that the outer edge of the inner disk is located at < 
0.21 AU and the inner wall of the outer disk is located at 30 AU from the central star (Espaillat 
et al. 2010, 2011). The disk around UX Tau A was spatially resolved by Sub-millimeter Array 
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at 880 fim wavelength, with the spatial resolution of 0.3" (Andrews et al. 2011). They found 
a dust-depleted disk cavity around the central star, and estimated the inner edge of the outer 
disk to be located 25 AU from the central star. UX Tau A shows no 10 /iin silicate emission 
(Espaillat et al. 2010), implying a lack of small dust grains in the disk. These features suggest 
dust grain growth in the circumstellar disk of UX Tau A. We conducted polarization imaging 
observations of UX Tau A in the if-band (1.6 /J,m) and investigate the collisional coalescence 
process of the dust grains in its disk. 

2. Observations 

Near-infrared (if -band; 1.6 fim) polarimetric imaging observations of UX Tau A were 
carried out 2009 December 23 with HiCIAO (High Contrast Instrument for the Subaru next gen- 
eration Adaptive Optics; Tamura et al. 2006) and the adaptive optics system, A0188 (Hayano 
et al. 2010), mounted on the Nasmyth platform of the Subaru Telescope (Table 1). These 
observations were conducted as part of the larger SEEDS survey (Tamura 2009). We em- 
ployed the Polarization Differential Imaging (PDI) mode. The Wollaston prism installed in 
HiCIAO divides incident light into two linearly polarized components which are perpendicular 
to each other and which are imaged simultaneously on the detector. Each image has 1024 x 
2048 pixels with a field of view of 9.75" x 20.09" and pixel scale of 9.521 mas/pixel in the 
east-west direction and 9.811 mas/pixel in the north-south direction. In the PDI mode, when 
the half-wave plate is set at the offset angle of 0°, 45°, 22.5°, and 67.5°, we obtain polarimetric 
images with the polarization direction at 0° and 90°, 90° and 0°, 45° and 135°, and 135° and 
45° components respectively. We used a coronagraphic mask with 0.3" diameter to suppress the 
brightness of UX Tau A. We obtained 44 frames, i.e. 11 frames per half-wave plate position, 
of UX Tau A with an exposure time of 60 s. The outside of the coronagraphic mask was not 
saturated. During the observations we fixed the star at the center of the coronagraphic mask, 
i.e. dithering method was not employed. In these frames, the tertiary component (UX Tau 
C) was also imaged outside the coronagraphic mask. The PSF reference star SAO93770 was 
imaged before observing UX Tau, using the PDI mode with the 0.3" diameter coronagraphic 
mask. The natural seeing was 0.5 - 0.6" in the i^-band. The A0188 measures the wavefront 
distortion by the atmospheric turbulence in the i?-band wavelength and compensates it at all 
wavelength. The FWHM of UX Tau C was 10 pixels (~ 0.1") on average. In the observations 
of SAO93770, with the ND filter in the AO device, we set the brightness difference between it 
and UX Tau to be 0.07 mag in the i?-band in order to level the AO correction performance 
in the observations of the object and the reference star. We took 11 frames with an exposure 
time of 90 s without employing dithering method. We obtained 10 flat frames on the same day 
as the target observations and 49 dark frames in 2009 December 24. Each exposure time was 
80 s for the flat frames and 5.6 s for the dark frames. The dark frames were used only for the 
identification of hot pixels. 
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Table 1. Properties of the observing objects 



Target RA [J2000] DEC [J2000] Sp.Type* H [mag]** Airmass 

UX Tau A 04 30 03.991 +18 13 49.39 K2 8.0 1.2 - 1.6 

UX Tau C 04 30 03.991 +18 13 49.39 M5 10.9 1.2 - 1.6 

SAO93770 04 07 08.734 +10 47 58.86 F8 7.0 1.2 - 1.3 

* Kraus & Hillenbrand (2009) for UX Tau A and (J. The Hipparcos and Tycho Catalogues for 
SAO93770. ** Corrcia et al. (2006) for UX Tau A and C. NOMAD catalog for SAO93770. 



3. Data Reduction 

The object frames were calibrated with the Image Reduction and Analysis Facility 
(IRAF). All frames of HiCIAO/A0188 have artifacts of horizontal stripes and vertical band- 
ings. The horizontal stripe has a size of 2048 x 64 pixels. Each frame has 32 horizontal stripes. 
We attribute the origin of these artifacts to a bias instability caused by temperature fluctua- 
tions between the detector and the Application Specific Integrated Circuit (ASIC). The ASIC 
translates the detector's analog signals into the digital signals. The vertical bandings occur at 
one pixel interval. They show an alternate direction between even and odd horizontal stripes. 

To mitigate these artifacts, first, the horizontal stripes were removed from the dark 
frames. We measured a median count for each horizontal stripe and subtracted it from each 
horizontal stripe. For the vertical bandings, even horizontal stripes were vertically flipped. We 
median-combined 32 horizontal stripes to make a master horizontal stripe. We subtracted the 
master stripe from each odd horizontal stripe. From each even horizontal stripe, the vertically 
flipped master stripe was subtracted. Finally, the dark frames were median-combined. This 
combined dark frame was only used as a mask frame for the extraction of hot pixels. 

Next, we eliminated hot and bad pixels from the object and flat frames. We regarded 
the pixels with more than 100 ADU in the combined dark frame as hot pixels. It is nearly 
equivalent to 12<t above the median of the combined dark frame. Hot pixels were replaced with 
10000 ADU and the other pixels were replaced with ADU. Using this frame as a hot pixel 
mask, hot pixels were interpolated by nearest good pixels in the object and flat frames. Then, 
we removed the stripe patterns from the object frames with the same procedure we used for 
the dark frames. In this process, however, we put circular masks on the stars and computed 
the median values of the stripes and the bandings except the mask regions. We considered that 
dark counts were also subtracted from the object frames through this destriping process. We 
did not subtract the stripes from the flat frames. Because the flat frames had counts two orders 
of magnitude larger than the stripes, we ignored the stripe modulations. The flat frames were 
median-combined and normalized to 1 ADU. The object frames were divided by the flat frame. 
Bad pixels and cosmic rays were interpolated by nearest good pixels in the object frames. 

After those processes, we obtained the flux images of the polarimetric components, F o 
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and -F 90 °, F go ° and F o, F 45 o and -F 135 o, and -F 135 o and -F 45 o, separated into left and right of 
the image by ~ 1070 pixels with each other. We cut the object frames into the left and right 
images. Each image has a 9.75" x 20.09" field of view. For each image, distortions caused by 
the Wollaston prism was corrected. Then, the central coordinates of the companion, UX Tau 
C, in all object images were measured. Using those coordinates, the object images were aligned. 
We also carried out the same procedures for the PSF reference star frames. The position of 
the photo-center of the reference star was estimated by fitting the halo component with a 2D 
Gaussian profile. Two frames of UX Tau A and three frames of the PSF reference star were 
not included in the following process, since they had insufficient tip-tilt correction. 

Linear polarization is described by the degree of polarization P and polarization angle 
a. These values are derived with the Stokes parameters, /, Q, and U 



90° 



I Fqo + Fg o 

U F 45 o — F 135 o 



/ -F45 + -^135° 



a 




(1) 
(2) 

(3) 
(4) 



where Fqo, i<9o°, -F450, and -F1350 are the intensity of each polarimetric component. By subtract- 
ing the right images from the left images, we derived 11 images of Q and —Q, and 10 images 
of U and —U. This subtraction method suppressed speckle noises that appeared both in the 
right and left images. We subtracted the — Q images from the Q images and the —U from the 
U. Thus, different aberrations between the left and right images due to the different dividing 
directions of the Wollaston prism were cancelled (double- difference technique; Hinkley et al. 
2009). 2Q and 2U images were divided by 2, then the Q and U images were combined respec- 
tively. Then, we corrected the instrumental polarization induced from the telescope and the 
instruments. We refer to Joos et al. (2008) for the correction. First, we construct a 2 x 2 Jones 
matrix for each optical element. In the HiCIAO case, we consider that the tertiary mirror and 
three mirrors of the derotator are probably dominant sources of the instrumental polarization. 
The surface materials of the mirrors expressed by the complex reflection index (n = 77 + in) 
was actually measured by the ellipsometer. We also describe the rotation of Stokes vectors as a 
function of the parallactic angle. Then, the multiplication of all single Jones matrices leads to 
a Jones matrix of whole telescope/instrument for a given set-up. Finally, the constructed Jones 
matrix was inverted and multiplied with the measured Stokes vectors to obtain the original 
Stokes vectors on the sky. In above method, we assumed monochromatic light and infinite 
F-number. 
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These images had distortions by the compensator in HiCIAO, the AO system, and 
the telescope. The compensator is located in front of the focal mask and reduces chromatic 
aberrations. By using images of M15 observed by the HST and HiCIAO, we calibrated the 
image distortion of Q and U. Multiplying both sides of equation (3) by /, we obtain the 
following equation, 



We constructed polarized intensity (PI) images with the Q and U images. Assuming that 
polarization due to the interstellar materials in front of UX Tau is negligible (See Section 4), 
the polarized intensity image represents the polarized component of circumstellar structures 
around the central star. 

In order to derive the polarization P of the circumstellar structures, the PI image 
needs to be divided by the I image which contains only the components of the circumstellar 
structures. The I image of UX Tau A (hereafter 7 tot ) consists of the intensity of the circumstellar 
structures as well as that of the central star. By subtracting the / component of the central 
star (hereafter /*) from the I to t, we obtain the / image of the circumstellar structures (hereafter 
-fdisk), i-e. /disk = hot — I*- We adopted the I image of the PSF reference star (hereafter J psf ) 
as the P. First, we constructed I images. We used the left and right images, F o and Fqq°, 
i<9o° and Fqo , F450 and -F1350 , and -F1350 and F450 . For UX Tau, we combined the left and right 
images into 7 tot images. For the PSF reference star, the central star position was measured in 
each image with fitting the halo component with a 2D Gaussian profile, then the position offset 
was corrected. Then we combined the left and right images into the I ps f image. The image 
distortion was also calibrated. Two frames of the PSF reference star have distorted PSFs, so 
that we used six frames of the PSF reference star in the following procedure. 

Since the Subaru Telescope is an alt-azimuth telescope and HiCIAO is mounted on the 
Nasmyth platform, the position angle of the spider changes in each image. The 7 ps f images 
were rotated to adjust their spider directions to those of every J tot image. The six J psf images 
of the PSF reference star were median-combined to make a PSF template / image (hereafter 
7*tmp) for each J tot image. Finally, we made 42 /*tm P images. From these J tot and i* tmp images, 
sky levels were subtracted. The sky value was the average of the mean counts of four regions. 
The regions were selected as the 100 x 100 pixel regions separated at least 5" from the central 
star. Next, we made the Jdi S k image. In each J tot and i* tmp image, we measured halo intensities 
of the central star. The evaluated regions were four regions of 20 x 20 pixels, which had a 
distance of 150 pixel (1.47") from the central star and did not overlap the spider directions. 
We used the ratio of the mean values of these halos to determine the proper scaling to apply 
to J*tm P image during PSF subtraction. We also shifted -Z* tmp image slightly so that intensity 
of the circumstellar disk would not show strong asymmetry in the residual image. We assume 
that the disk intensity is symmetric, thus we may have over-corrected the PSF subtraction. 




(5) 
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The amount of the shift was 1.3 pixels (0.01" = 0.1 FWHM ~ 2 AU) on average. Therefore, 
the inner working angle of the /disk image is ~ 23 AU. Then, we subtracted /* tmp image from 
/tot image. This process was performed in each / tot image. As a result, we obtained 42 /disk 
images containing only the circumstellar structure. They were median-combined. Finally, the 
PI image was divided by the /disk image, thus the P image of the circumstellar structure was 
generated. 

4. Results 

We detected strongly polarized components around UX Tau A. The polarized intensity 
(PI) image shows an elliptical structure (Figure 1). We evaluate the extent of the region where 
the intensity is more than five times the standard deviation in the sky region of the PI image. 
It has a semi-major axis of ~0.83" (~ 120AU) with the position angle of 165° ± 2° and a 
semi-minor axis of 0.58" (~ 80AU). 

We associate this polarization structure to the circumstellar disk of UX Tau A. If the 
circumstellar disk of UX Tau A has a circular and geometrically thin structure, the disk is tilted 
to east- west with the inclination of 46° ± 2°. The west side of the structure is brighter than the 
east side. We assume that the west side of the structure is the near side to us. We measured 
degrees and angles of polarization in every 9x9 pixel region. The degree of polarization in 
regions with SNR > 5<r ranges from 1.6 to 66 %. The polarization vectors show a mostly 
centrosymmetric pattern centered on UX Tau A, supporting our claim that this structure is 
attributable to its circumstellar disk. The radial profiles of the polarized intensity along the 
semi- major axes are shown in Figure 2. The polarized intensity declines with the third power 
of the disk radius. 

Figure 3 shows an azimuthal profile of the polarization degree (P) along the ellipse of 
e = 0.72 (i = 46°) and a semi-major axis of 0.30" (42 AU). We measured mean polarization 
degrees and standard deviations in every 9x9 pixel region. It is sinusoidal; two maxima and 
two minima appear with 90° intervals. This pattern is prominent up to the semi-major axis 
of ~ 0.5". Figure 4 represents radial profiles of the polarization degree along the semi-major 
axes. We measured mean polarization degrees and standard deviations in every 9x9 pixel 
region. The degree of polarization declines with radial distance to ~ 0.5". Beyond 0.8" from the 
central star, the polarization is not significantly detected. Assuming that the polarized intensity 
and the polarization degree of the disk are constant within binned 9x9 pixels respectively, 
these uncertainties are up to ~ 80 ADU and ~ 20 % as shown in Figures 2 and 3. These 
estimates do not take into account any systematic effects from the data reduction such as the 
PSF subtraction and the correction for the instrumental polarization. Further discussion on 
the instrumentation and the systematic errors will be presented elsewhere. 

Polarization due to interstellar medium in front of UX Tau binary system is negligible. 
The system is slightly embedded in a molecular cloud; the interstellar extinction of UX Tau A 
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Fig. 1. A polarized intensity image of UX Tau A. The polarization vectors are superimposed in the 
region where the polarized intensity is larger than 5o\ The vector lengths and directions indicate the 
degree and angle of the polarization, respectively. The vector scale is shown in the bottom- left. A circle 
shows the 0.3" diameter coronagraphic mask positioned on the central star. The bottom color-bar denotes 
the polarized intensity in ADU. 

is 1.8 mag in the ^-band (Espaillat et al. 2010). Serkowski et al. (1975) obtained a relationship 
between polarization and reddening as E(B — V) ^ P/9. We would therefore anticipate little 
foreground interstellar polarization would be present. Indeed, the polarization degree of UX 
Tau C, which was measured in 50 x 50 pixels centered on the star, is only 0.5 %. We do not 
detect any circumstellar structures around UX Tau C. Since no 880 /mi emission is detected 
from UX Tau C (Andrews et al. 2011), this source is unlikely to have circumstellar structures 
and its measured polarization can be assumed to be a reasonable proxy for the interstellar 
polarization along the line of sight to UX Tau A. Therefore, the interstellar medium in front of 
the UX Tau system has only negligible polarization. 

Outflows emanating to the north-west and south-east directions are also possible sources 
to explain the elliptical structure. Lucas et al. (2004) detected the polarization structure 
around HL Tau and suggested that its strong polarization structure is caused by bipolar out- 
flows. However, the disk around UX Tau A has already been imaged by Sub-millimeter Array 
(Andrews et al. 2011). The inclination and the position angle of the semi-major axis of the 
disk were determined to be 35° and 176°, respectively. The disk geometry is roughly consistent 
with that derived from the polarized intensity image. Therefore, we consider that the elliptical 
structure we detected is not attributed to the outflow, but to the circumstellar disk around the 
central star. 
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Fig. 2. Radial profiles of the polarized intensity measured along semi-major axes of the circumstellar 
disk of UX Tau A. The plots and the errorbars show mean polarized intensity and standard deviation in 
every 9x9 pixel region. The filled and open circles denote the polarization intensity along semi-major 
axes in the north-west and south-east directions, respectively. The solid line show a power law fitting to 
the observed radial profile. The polarized intensity declines as the third power of the disk radius. 




Fig. 3. An azimuthal profile of the polarization degree of the circumstellar disk. It is measured on the 
ellipse of eccentricity e = 0.72 and a semi-major axis of 0.30" and 0.50", which corresponds to the disk 
radius of ~ 40 AU and 70 AU. For the position angle (PA), north is 0°and east is 90°. The plots show 
mean polarization degrees binned in 9 x 9 pixels, and maxima and minima are clearly seen in 90° intervals. 
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Fig. 4. Radial profiles of the polarization degrees measured along semi-major axes of the circumstellar 
disk of UX Tau A. The plots show mean polarization degrees binned in 9 x 9 pixels. The filled and open 
circles denote the polarization degrees along semi-major axes in the north-west and south-east directions, 
respectively. From 0.15" to 0.50", the degrees of polarization show decreasing trends. 

5. Discussion 

5.1. Geometry of the circumstellar disk 

The circumstellar disk looks smooth. A gap-like structure at 0.21 - 30 AU is suggested 
from the disk model fitting to the SED (Espaillat et al. 2010, 2011). A disk cavity is also 
resolved by the SMA observations (Andrews et al. 2011), which is identical with the outer 
radius of the gap structure. However, we do not find such a structure. Our data do not show 
any signs for a gap or cavity down to the achieved inner working angle (~ 23 AU). We are not 
sensitive to a possible inner disk with the radius of < 0.21 AU (Espaillat et al. 2010) because 
it is within the coronagraphic mask. 

We examine whether UX Tau C influences the disk geometry of UX Tau A. The 
Lagrangian point LI is the point where gravitational forces of the primary and tertiary stars 
balance each other to a mass point. With the primary mass of 1.3M , the tertiary mass of 
O.16M (Kraus & Hillenbrand 2009), and the separation of 2.63" (White & Ghez 2001), the 
location of the point LI was derived to be 1.9" from the primary star. Given that the disk of 
UX Tau A extends up to 0.8", we conclude that the disk is located in the gravitational zone of 
the primary star and is not gravitationally affected by UX Tau C. 

5.2. Polarization models of the circumstellar disk 

We construct several polarization models of the circumstellar disk. First, we define 
the disk geometry. A scattering angle is an angle of deflection from the forward direction of 
incident light. We consider the disk model based on McCabe et al. (2002). In the disk geometry 
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of Figure 5, the scattering angle is calculated as 




1 + cos 2 (PA -PA)) tan 2 i 



1 




(6) 



where i is the disk inclination to the line of sight, PA is the position angle measured on the 
mid-plane of the disk, and PAq is the position angle on the near side of the semi-minor axis 
of the disk. ope n is the opening angle of the disk. It expresses the disk height with optical 
thickness r > 1. If cosPA > 0, j = 1, and if cosPA < 0, j = (McCabe et al. 2002). Substituting 
the inclination, the semi-minor axis, and the height of the disk around UX Tau A into equation 
(6) allowed us to acquire the scattering angle at each position angle. The disk inclination of 
46° was applied. From the semi-major axis of 165° (Figure 1), we assumed that the semi-minor 
axis direction was 255°. 

Next, we consider three scattering approximations. The size parameter for dust grains 
is defined as X = 1i\aj\. X = 1 corresponds to 0.25 /xm of the grain radius (a) for the P-band 
observations. We first adopt, as the simplest case, a Rayleigh scattering approximation. The 
Rayleigh scattering approximation holds in the case of sufficiently small dust grains compared 
to the observing wavelength (X <1). For the P-band observations, it can be applied for the 
grains with the radii (a) much smaller than 0.25 /im. The relationship between the degree of 
polarization and the scattering angle is represented by the following equation (van de Hulst 



The degree of polarization P shows a maximum value P max at 6 scat = 90° and a minimum value 
Pnin at scat = 0° and 180°. By combining equations (6) and (7), we are able to calculate the 
polarization pattern on a circumstellar disk. 

We consider a geometrically thin disk (0 opC n = 0°). It produces a sinusoidal azimuthal 
profile of the polarization degree (Figure 6). Maximum polarization degrees are ~ 60 % at 
PA(P max ) ~ 165°, 345°. Minimum degrees are - 20% at PA(P min ) ~ 75°, 255°. Two P max 
and two P m i n have the same values respectively. The interval between PA(P max ) and PA(P mm ) 
is 90°. Although the data points match the observational polarization degrees quite well for 
P max , they shows large differences (about 20 %) from the observational degrees for P mm - The 
case of a geometrically thick disk with ope n = 30° is also shown in Figure 6. Whereas two 
Pmax have the same values, two P m m have different values. Thus, we conclude that neither disk 
model matches the observational azimuthal profile of the polarization degrees. Espaillat et al. 
(2011) suggested that maximum grain size is 10 /im at the inner edge of the inner disk and 5 
jum at the inner edge of the outer disk, from fitting a disk model to the SED of UX Tau A. If 
the disk consists only of such large dust grains, the Rayleigh scattering approximation is not 
appropriate for the P-band observations. 

Next, we consider a Mie scattering approximation, which is applicable for spherical dust 



1957) 



P 



1 - cos 2 6> scat 
1 + cos 2 6> scat 



(7) 
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grains with the size comparable to the observing wavelength (X « 1, i.e. a ~ 0.25 /mi). The 
degree of polarization depends not only on the scattering angle, but also on the grain size, 
the size distribution, the refraction index, and other grain parameters. Using the observing 
wavelength, the grain size, the grain size distribution, the refraction index (n), and the absorbing 
coefficient (k), the Mie scattering program allows us to obtain the degree of polarization as a 
function of the scattering angle (Figure 7). A similar work done on comet theory using Mie 
theory, details the procedure (Sen et al. 1991). We found that large dust grains show low P max . 
In the case that the largest dust grains have a radius as large as 1 mm, the polarization degree 
shows negative for the large scattering angle (Murakawa 2010). The negative polarization 
indicates a radial pattern of the polarization vectors. We conclude that any Mie scattering 
approximation does not reproduce the observed polarizations of UX Tau A. Even if we use a 
geometrically thick disk, the models are hardly different from that of a geometrically thin disk. 

Finally, we focused on geometric optics. It is applicable for interpreting light scattering 
by irregular shaped particles sufficiently larger than the observing wavelength (X 3> 1, i.e. a ^> 
0.25 /mi). We adopted computer simulation results of Grynko & Shkuratov (2003), in which 
light scattering is calculated for various grains. We used the polarization profiles for 100 faceted 
particles. Grynko & Shkuratov (2003) calculated it with the refractive indices of some silicate 
(n = 1.5, k = 0.004), and a ~ 30 /mi. We considered that the 100 faceted particles presumably 
correspond not to the dust grains with smooth spheres but to those with rugged surfaces. The 
degree of polarization shows % at scat = 0°, 180° and maximum value at 6> scat ~ 85°. These 
distribution of the polarization degree resembles that of the Rayleigh scattering approximation. 
One characteristic is that the degree of polarization falls down to about 10 % at # scat ~ 45° and 
135°, whereas that in the Rayleigh scattering approximation has about 30 % at # scat = 45° and 
135°. With this distribution of the polarization degrees, we constructed the polarization model 
of a geometrically thin disk (Figure 8). Maximum polarization degrees are ~ 80% at Pv4(P max ) 
= 165° and 335°. Minimum degrees are ~ 3 % at PA(P min ) = 55° - 85° and 225° - 275°. 
Although the model shows maximum degrees about 10 % larger than the observed maximum 
degree at PA ~ 165°, it reproduces well the observed minimum polarization degree. Moreover, 
PA(P max ) and PA(P min ) of the model are consistent with those of the observed PAs. We did 
not obtain such a profile in the Mie scattering approximation by spherical grains with a = 30 
//hi. Therefore, we suggest that dust particles in the disk are nonspherical grain with 60 /mi 
diameter at 40 AU from the central star. This argument is consistent with no 10 /im silicate 
emission in the UX Tau A spectra, indicating a lack of small dust grains (Espaillat et al. 2010). 

We also considered a geometrically thick disk. Assuming vertical hydrostatic equilibrium 
and no turbulence in the disk, the density profile of the disk is given by 



where H is the scale height and po is the density at the disk mid-plane. Muto (2011) indicates 




(8) 



14 



that the disk thickness Z with r = 1 corresponds to 3H for such a disk. The large opening 
angle of the circumstellar disk, op cn, corresponds to large ratio of the scale height to the disk 
radius (H/R). We simulated thick disks with various O pen values. A disk with large opC n value 
shows large interval of PA(P max ). For example for the disk with H/R of 0.1, PA(P max ) results 
in 150° and 350°. Moreover, minimum polarization degrees at PA(P m [ n ) = 70° and 250° have 
different values. It is distinctly different from the polarization profile of UX Tau A. If <f) p Cn < 
11°, corresponding to H/R < 0.067, the models are consistent with the polarization profile of 
UX Tau A. We conclude that the geometrical thickness of the UX Tau A disk agrees with H/R 
= - 0.067. Most of the dust grains in the disk have therefore settled into the mid-plane and 
formed the flat disk. Andrews et al. (2011) found a disk geometry favoring a rather flat disk 
structure. 

The polarization degrees decrease with the semi- major axes (Figure 4). Since the 
Keplerian velocity is small in the outer part of the disk, it is generally considered that the 
dust particles in the outer disk have small relative velocities. We consider a slow growth rate 
of such dust, thus resulting in small dust grains in the outer part of the disk. However, the low 
polarization degrees in the outer part of the disk are not reproduced by the Rayleigh scattering 
approximation. The computer simulations of light scattering with geometric optics approxi- 
mation (Grynko & Shkuratov 2003) calculated the polarization degree profile for grains with 
various number of facets. They indicated that the maximum polarization degree decreases with 
decreasing the facet number. We expect dust grains with rough surface in the outer part of 
the circumstellar disk. Most previous studies of scattered light images from circumstellar disks 
could explain the scattering (or polarization) function with small, ISM-like grains. For AB 
Aur, for instance, Perrin et al. (2009) found a similar P max as in this present study, but they 
could explain the results with more or less typical, small ISM grains. Only recently, Quanz 
et al. (2011) found some indications based on PDI data, that on the surface of the HD 100546 
disk micron sized grains might be present. Andrews et al. (2011) used also a rather ISM-like 
dust population to explain both their interferometric observations and the SED of UX Tau A. 
However very small P m ; n observed in this work can be reproduced by the scattering model of 
the dust grains with rough surfaces. To obtain further information of shapes and sizes of the 
grains, geometric optics models are to be investigated with various parameters. 

5.3. Dust growth in the circumstellar disk 

The forces acting on dust grains in a protoplanetary disk are gas drag force and grav- 
itational force of the central star. In a condition with a temperature typically expected in a 
circumstellar disk, grains with the sizes less than a few tens micrometers are marginally affected 
by gas drag forces and dominated by Brownian thermal motion (Weidenschilling 2000). As a 
consequence of frequent collision and coagulation, the dust grains grow larger. We examine 
how large dusts can grow during settling into the mid-plane at 40 AU from the central star. 
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It is claimed that the dust grains grow to several centimeters at 1 AU, during sedimentation 
(Nakagawa et al. 1981). We calculated the radii of grains settling toward the mid-plane accord- 
ing to the equations for dust growth during sedimentation (Takeuchi et al. 2009). We assumed 
the dust sticking probability C = 1 and the grain bulk density of 1 g cm -3 . We do not consider 
any turbulence in the disk. As a result, we found that the dust grains can grow up to 100 /mi 
radius during sedimentation at 40 AU from the central star. Therefore, it is possible for dust 
grains to grow to 30 /mi radius at 40 AU from the central star. 

Next we evaluate the timescale that the dust grains grow from 0.1 /mi to 100 /xm during 
sedimentation toward the disk plane. We take account of the gas drag. If the circumstellar disk 
of UX Tau A is optically thin at sub-millimeter wavelengths, the sub-millimeter observations 
(Andrews et al. 2011) indicate that the gas surface density of the UX Tau A disk at 40 AU from 
the central star is ~ 10 g cm -2 . For the disk with H/R = 0.067, the disk thickness with r = 1 
is about 16 AU at the disk radius of 40 AU. We obtain the average gas density of 4.2 x 10~ 14 
g cm -3 (2.5 x 10 10 atom cm -3 ). With the relationship between the gas density and the mean 
free path of gas molecules (Takeuchi et al. 2009), we derive 4.8 x 10 4 cm for the gas mean free 
path. Since the radii of the dust grains up to 100 /mi is sufficiently smaller than the gas mean 
free path, the dust grains grow under the Epstein regime. As a consequence of calculating the 
equation for the sedimentation timescale (Takeuchi et al. 2009), it is revealed that it takes at 
least 10 5 years for dust grains with 0.1 /mi radii to settle and grow up to 100 /mi in radius. 
The age of UX Tau A is estimated to be about 10 6 years (Kraus & Hillenbrand 2009). We 
claim that the coagulated growing dusts have already settled into the equatorial plane of the 
circumstellar disk. It is consistent with our observational evidence that the circumstellar disk 
of UX Tau A is geometrically thin. 

Dust grains with sizes between a few micrometer and a hundred micrometer are entrained 
by gas toward the central star. Adachi et al. (1976) investigated the spiral motion of the dust 
grains in a protoplanetary disk. The decay time of the spiral motion corresponds to the infall 
timescale of the dust grains toward the central star. Applying their calculations, with the 
radially constant gas density of 4.2 x 10~ 14 g cm" 3 , we derived 2 x 10 6 years for the dust 
grains with radii of 100 /mi to infall from 40 AU to the central star. We therefore suggest that 
the circumstellar disk of UX Tau A still contains such large dust grains. 

6. Conclusions 

With HiCIAO/A0188 mounted on the Subaru Telescope, we carried out if -band po- 
larimetric imaging observations of UX Tau A, which has been classified as a pre-transitional 
disk object. The observation revealed a circumstellar disk around UX Tau A at the spatial 
resolution of 0.1" beyond 23 AU from the central star. The disk extends to 120 AU in radius. 
The disk is inclined by 46° ± 2° as the west side is near. We did not detect evidences of the 
inner disk and the gap-like structure. 
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Fig. 5. Schematic view of the disk geometry. 
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Fig. 6. Azimuthal profiles of the polarization degree. The filled circles show the observed degrees of 
polarization measured along the ellipse with an eccentricity of 0.72 and a semi-major axis of 0.30". The 
solid line represents the polarization degree of a geometrically thin disk model (<^> pen = 0°) and the dashed 
line shows that of a geometrically thick disk model of (f> op en = 30°. We used the Rayleigh scattering 
approximation. The polarization profiles of the models are convolved with a 20° window in order to match 
the spatial resolution of the 9x9 pixels binned profile of the observational polarization. 
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Fig. 7. Azimuthal profiles of the polarization degree. The filled circles show the observed degrees of 
polarization measured along the ellipse with an eccentricity of 0.72 and a semi-major axis of 0.30". We 
constructed four geometrically thin disk models with the maximum grain radius of 0.25, 0.50, 1.0, and 
10 /im, using the Mic scattering approximation. The polarization profile of each model is shown by the 
dot, the short dashed, the long dashed, the dashed dotted lines, respectively. We used 0.005 /jm for the 
minimum grain radius and N(a) oc a~ 35 for the grain size distribution. The polarization profiles of the 
models are convolved with a 20° window in order to match the spatial resolution of the 9x9 pixels binned 
profile of the observational polarization. 



model open ^O° 

model opea = 1 1 ' 

• observed 




50 100 150 200 250 300 350 



PA [deg 



Fig. 8. Azimuthal profiles of the polarization degrees of the circumstellar disk. The filled circles show the 
observed degrees of polarization measured along the ellipse with an eccentricity of 0.72 and a semi- major 
axis of 0.30". We showed two disk models constructed by the geometric optics approximation. The solid 
line represents the geometrically thin disk model. The dashed line shows the geometrically thick disk 
model with 4> opcn = 11°. The polarization profiles of the models are convolved with a 20 window. Both 
models approximately reproduce the observed polarization degrees. 
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The notable feature of the circumstellar disk is the huge variation of the polarization 
degrees. It varies from 1.6 to 66 %. We constructed several polarization models of the cir- 
cumstellar disk based on the Rayleigh scattering and Mie scattering approximations. However 
neither models were consistent with the observational azimuthal profile of the polarization 
degrees. Focusing on geometric optics, we built the polarization model of the geometrically 
thin disk with nonspherical grains with the radii of 30 jum. The model reproduced well the 
observational azimuthal profile of the polarization degree. We suggest that UX Tau A has a 
geometrically thin disk containing the nonspherical dust grains with the radii of 30 /xm. 

Such a disk with nonspherical large dust grains is consistent with the core accretion 
model of planetary formation process. At 40 AU from UX Tau A, the dust grains can grow 
up to 100 fim in radius by collisional coagulation and settle toward the mid-plane with the 
timescale of 10 5 years at the earliest. Observational evidence of large dust grains as well as 
the gap structure in the circumstellar disk provides robust signatures of planetary formation 
process in the UX Tau A system. 
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